Supersolid phases in the one dimensional extended soft core Bosonic Hubbard model 
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We present results of Quantum Monte Carlo simulations for the soft core extended bosonic Hub- 
bard model in one dimension exhibiting the presence of supersolid phases similar to those recently 
found in two dimensions. We find that in one and two dimensions, the insulator-supersolid transi- 
tion has dynamic critical exponent z = 2 whereas the first order insulator-superfluid transition in 
two dimensions is replaced by a continuous transition with z = 1 in one dimension. We present evi- 
dence that this transition is in the Kosterlitz-Thouless universality class and discuss the mechanism 
behind this difference. The simultaneous presence of two types of quasi long range order results in 
two soliton-like dips in the excitation spectrum. 

PACS numbers: 75.40.Gb, 75.40.Mg, 75.10.Jm, 75.30.Ds 



Experimental and theoretical work on supersolidity 
(SS), a thermo dynamically stable phase in which diag- 
onal (density) and off-diagonal (superfluid) order coex- 
ist Q, has seen an enormous resurgence over the past 
several years. In the case of solid Helium, the obser- 
vation of an anomalous moment of inertia |2| has been 
vigorously re-investigated by several groups both experi- 
mentally and theoretically . In parallel, an entirely 
new arena for the investigation of supersolid phenomena 
has arisen in the context of cold atoms in optical lattices, 
whose description by the bosonic Hubbard Hamiltonian 
is well established |5(|. It is known that this model does 
not exhibit supersolid phases due to the absence of longer 
range interactions. However, the recent realization of 
Bose-Einstein condensate of 52 Cr atoms, which have an 
exceptionally large magnetic dipole moment and there- 
fore long-range interactions, has opened wide the possi- 
bility of experimantal observation of SS phases on con- 
fined optical lattices Q. 

Such strongly dipolar ultra-cold atoms, when placed 
on optical lattices, are governed by the extended Bosonic 
Hubbard Hamiltonian which has been shown, via exten- 
sive Quantum Monte Carlo (QMC) simulations, to ex- 
hibit SS phases in two dimensions. The model we con- 
sider is, 



H = -t (at a 3 + a j a i) + -2^2 n i( n i - !) 

<ij> i 

+V n i n 3- (!) 

<ij> 

The hopping parameter, t, sets the energy scale, = 
a\a,i and [a,*,aj] = <5y are bosonic creation and destruc- 
tion operators. The repulsive contact (near neighbor) in- 
teraction is U (V). Recently, it was shown || in 2d that 
the model described by Eq. has a SS phase if U and 
V are such that when the charge density wave (CDW) 
phase is doped above half filling, the added bosons go 
on already occupied sites (roughly U < 4V). Doping 
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FIG. 1: The boson density, p, vs the chemical potential, fx, 
for L — 64, j3 — 20 with 200 time slices. Inset: p 3 vs p as 
p — 1/2 is approached from below. See text near Eq. ©. 



below half filling provokes phase separation and a first 
order transition to a SF phase. If multiple occupancy is 
suppressed, there is phase separation when the system is 
doped above half filling as in the hard core case 9|. This 
raises the possibility that an analogous phase exists for 
the one dimensional form of Eq. Q in which case this 
would be the first example of a SS in a one dimensional 
system. 

In this letter we re-examine the one dimensional soft 
core bosonic Hubbard model 

H3E2, Eq. © under the 
same conditions. We conclude that a ground state SS 
phase does exist with similarities to, and important dif- 
ferences from, that in 2d. Calculations of the excitation 
spectra and Landau critical velocity demonstrate that 
the SS phase is stable, and allow us to infer the dynamic 
critical exponents. We use the World Line algorithm in 
our QMC simulations. 

We begin by demonstrating that for V sufficiently 



2 



large, the system has gapped, incompressible CDW in- 
sulating behavior at half-filling, p = 1/2. Figure ^ shows 
p vs fi for U = Wt and V = 6t, 7t, 8t, 9t using the WL 
algorithm. The plateau at p = 1/2 clearly indicates a 
non-zero gap for all four cases. However, the behavior of 
the system depends on whether it is doped below or above 
p = 1/2, and on the value of V. Starting with the CDW 
and removing a particle creates a "defect" with three ad- 
jacent unoccupied sites which then splits into two mobile 
two-hole bound states (solitons) which destroy the long 
range CDW order in favor of quasi-long range order with 
p s 7^ 0. Consequently, all the curves coincide for p < 1/2 
in Fig^ In contradistinction, when the 2d CDW phase 
is doped below half filling, phase separation takes place 
leading to a first order phase transition. Upon adding a 
particle to the CDW for V < 6t, the extra particle goes to 
an unoccupied site since the cost of multiple occupancy 
exceeds that of near neighbors. Thus, a defect of three 
adjacent occupied sites is created which then splits into 
two mobile two-particle bound states (solitons) which are 
the particle-hole symmetric analogs of the p < 1/2 case. 
This particle-hole symmetry, seen for V = 6t in Fig. ^ is 
of course approximate for finite U. On the other hand, 
upon adding a particle to the CDW when V > 6t, the cost 
of near neighbors exceeds that of multiple occupancy, and 
the extra particle goes to an already occupied site (which 
cannot happen in the hardcore case). This breakdown of 
particle- hole symmetry is seen for V = 7t, 8t, 9t in Fig.^ 
For p > 1/2 the curve first rises very sharply hinting at 
a diverging compressibility, k = dp/ dp which is not the 
case for p < 1/2. We shall return to this point below. 

In order to characterise the phases, we calculate the 
SF density, p s |15|. and the structure factor, S(k), 
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n(x)n(x')), 



(2) 
(3) 



where W is the winding number, L the length of the 
system, d the dimension and (n(x)n(x')) is the density- 
density correlation function. Figure^a) shows p s versus 
p for U = lOt and V = 6t, 7t, 8t, 9t while H>) shows 
S(k = 7r) for the same values. For p < 0.5, S(n) is essen- 
tially zero, i.e. there is no CDW, while at the same time 
p s =/= 0. This is a SF phase. This behavior is qualita- 
tively different from that in two dimensions 0, E| where 
the system undergoes phase separation and a first order 
transition. Such phase separation cannot happen in one 
dimension, and indeed Fig. ^] does not exhibit any nega- 
tive compressibility reg ions [9j. Instead, solitonic excita- 
tions are produced [12j. 

For p > 0.5 and V > 6t, S(ir) takes non-zero val- 
ues indicating simultaneous co-existence of SF and long 
range CDW order, in other words a SS in one dimension. 
That this SS does not undergo phase separation is again 
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FIG. 2: (a) The superfluid density, p 3 , and (b) the structure 
factor, S(k), as functions of the number density p for four 
values of V. L = 64 sites, f3 = 20 and U = Wt. 



clear from Fig. ^ where |9j the compressibility, k is never 
negative. For V — 7t, finite size scaling shows that in 
the thermodynamic limit, p s maintains a non-zero value 
whereas the behavior of $(71") depends on p as L — > 00. In 
other words, as this system is doped above half filling, it 
first enters a SS phase, then a SF phase and finally back 
into a SS. The phase diagram in two dimensions (Fig. 
1 in Ref. 0) does not exhibit such re-entrant behavior. 
On the other hand, for V = 8t, 9t, the system remains in 
the SS phase for all p > 0.5, at least up to the highest 
density we studied, p = 1.1, except at p = 1 where it 
is a gapped CDW insulator with alternating empty and 
doubly occupied sites. 

The phase diagram for a fixed value of U in the 
{p/t, t/V) plane is obtained by sweeping p/t at con- 
stant t/V as shown in Fig.Q] This is shown in Fig. E3 for 
U = lOt. The phases and the nature of the transitions 
between them are explained in the caption. The details 
of the phase diagram above p = 1 have not been deter- 
mined beyond establishing the natures of the phases just 
above the CDW-II and Mott-I as indicated in the figure. 
When CDW-II (Mott-I) is doped above or below p = 1, 
the system becomes supersolid (superfluid). We remark 
that in figure E|(C/ = lOt) the CDW-II and Mott-I lobes 
have a very small contact region. This contact region 
increases as U is increased, i.e. the gap at the transition 
increases. On the other hand, if U is decreased, the two 
lobes separate and the direct transition between CDW-II 
and Mott-I is lost: a SF phase intervenes. 

Having shown that supersolid phases exist, we turn 
now to the excitation spectrum, fi(fc) which, along with 
p s , characterizes the superfluid phase. If Q(k) oc k for 
small fc, excitations are phonons and the superfluid is 
said to be stable according to the Landau criterion: the 
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FIG. 3: The phase diagram for U = lOi. CDW-I is alternat- 
ing empty and singly occupied sites, CDW-II is alternating 
empty and doubly occupied sites, Mott-I is the first Mott 
lobe with p = 1. Solid lines are second order transitions with 
z = 2; thick solid lines have z = 1 and appear to be in the 
KT universality class (see Eq. ©); the CDW-II to Mott-I 
transition is first order. Dashed lines connect schematically 
the numerical results at the smallest t/V to the exact values 
at t/V = 0. 



Landau critical velocity, v^, is finite and is given by the 
slope of the line passing through the origin and tangent to 
the roton (or, for d = 1, the soliton) minimum. To obtain 
fi(A;), one can calculate, via QMC, the space and imagi- 
nary time separated density-density correlation function 
which then gives the dynamic correlation function. With 
the help of a numerical Laplace transform, for example 
using the Maximum Entropy method, Q (k) can be ob- 
tained 16). This procedure is demanding numerically. 
However, it was shown that the rigorous Maxent 
procedure is not necessary since one obtains a good ap- 
proximation (upper bound) for Q (k) by using the /-sum 
rule leading to Feynman expression, 



ft(k) 



Ek 
LS(k)' 



(4) 



where, 



Ek = — (cos k 
L 



L 

E 



(4 



a i+ i + a] +1 aA |* ) 



(5) 

and l^o) is the ground state of the system. 

Figure0]shows fl(k) vs k for several values of the den- 
sity and interaction strength. The first three cases listed 
in the legend correspond to a fixed density p — 15/32, 
just below half filling in the SF phase. We note: (a) fi- 
nite size effects are negligible, (b) f2(fc) oc k as k — > 0, 
indicating phonon excitations and a stable SF, and (c) 
the soliton minimum is clearly visible at k = 2.95. As 
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FIG. 4: The dispersion relation, £l(k), vs the wave- vector, k, 
for several system sizes and boson numbers Nb. For all cases 

U = m. 



p — ► 1/2 from below, Q(k) remains linear for small k but 
the soliton minimum moves towards k = n and its energy 
approaches zero, pinching the fc-axis, and signaling the 
CDW instability. Therefore, the Landau critical veloc- 
ity, vj? — > although the SF remains stable, with a high 
phonon velocity, all the way to the transition. The lin- 
earity of O(fc) at the transition shows that the SF-CDW 
transition has a dynamic critical exponent z = 1. 

The next three cases listed in the legend, Fig. 0] cor- 
respond to p = 17/32, just inside the SS phase. Again 
finite size effects are under control. In addition, we see 
that f2(fc) softens as k — > and is no longer linear and 
that this effect increases with increasing V (the last case 
in Fig. 0| . Fitting the form fl(k) oc k z to the small k dis- 
persion curves and sizes up to L — 96 for p = 0.51 gives 
z ~ 1.9. This suggests a dynamical critical exponent 
z = 2 for the CDW-SS transition. The same behavior is 
observed for fi(fc) as the p = 1 CDW-II phase (alternat- 
ing doubly occupied and vacant sites) is approached from 
the SS phase above or below full filling, p — > l^. There- 
fore, the quantum phase transitions to the p — 1/2 CDW 
are different depending on whether they are approached 
from below (f2(fc — > 0) oc k and soliton minimum pinch- 
ing the fc-axis) or from above (f2(fc) softens to a quadratic 
form). However, the transitions to the CDW-II phase, 
p — > l^, are in the same universality class with z = 2. 
We also have performed simulations in two dimensions 
and verified that the CDW-SS transition has z = 2. 

Therefore, Q(k) yields for the SF-CDW transition z = 
1 while z = 2 for the SS-CDW transitions both at p = 1 /2 
and p = 1 (and presumably also for p — 3/2 etc). In 
addition, it was shown in Ref. [I3 . lr7T | that, when V = 0, 
p s scales as p s ~ |p c ~ p| z_1 with z — 2 where p c — 1,2, .. 
is the density at which the MI phase is established. Our 
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FIG. 5: Comparison of the structure factor, S(k) for the same 
doping above and below half filling. L = 64, /3 = 20, U = 
lOt, V = 8t. Inset: The corresponding dispersion, Q(fc). 

data show that for the SS-CDW p s ~ |p c — p| 95 again 
giving z « 2 assuming the same scaling law. On the 
other hand, for the SF-CDW transition z = 1 yielding 
Ps ~ |Pc — p| (Pc = 1/2 here) which means that for 
this transition p s does not vanish as a power law. In the 
absence of power law scaling, it is natural to attempt a 
Kosterlitz-Thouless (KT) scaling form, 

Ps = Ae- B ^\ p - p ^ (6) 

where A and B are nonuniversal constants. This is il- 
lustrated in the inset of Fig. ^ which shows p s vs p for 
L = 64, U = lOt, V = 8.5t. The solid line is a two 
parameter fit with Eq. ©, the dashed line is a two pa- 
rameter power law fit to the same points. The goodness 
of this fit suggest that this transition is in the KT uni- 
versality class. 

The dispersion curves in Fig.0]show clearly that when 
the system is doped below p = 1/2, Q(k) develops a soli- 
ton dip. Figure [S] shows the origin of this behavior. For 
Nf, = 19, fi(fc) shows a pronounced minimum, while for 
2V& = 45 it shows two milder ones (inset Fig. [5) . These 
excitation dips are caused by peaks in the structure fac- 
tor which point to the presence of quasi long range order 
with the corresponding ordering vectors, k*. It is in- 
teresting to note that for Nb = 45 (and in general for 
p > 1/2) the structure factor S(k) has two peaks (Fig. EJ) 
indicating the simultaneous presence of two types of quasi 
long range order. It is clear from the figure that one of 
the peaks occurs at exactly the same k* as the soliton 
peak for the same doping below half filling. This peak for 
p > 1/2 is caused by bosons occupying previously empty 
sites thus aquiring near neighbors even though the large 
value of V tries to suppress this. In other words, this 
is the vestige of the approximate particle-hole symmetry 



present for U = lOt and V < 6t. The peak at the lower 
k comes from the quasi long range order produced when 
doubly occupied sites are produced. Recall that at full 
filling the ground state is a CDW with alternating doubly 
occupied and empty sites. 

We now return to the behavior of the compressibility, 
K = dp/ dp, as p — ► l/2 ± . It was argued in [lj[ for the 
SF-MI transition (V = 0, p = 1) that k ~ \p- p c \ vi ~ d ~ z \ 
where p c is the chemical potential at the transition. For 
the p — > l/2 + transition, {d = 1, z — 2), our preliminary 
results indicate k ~ |/x — p c \~°' 6 which is consistent with 
v = 1/2 assuming the same scaling. For the p — > l/2~ 
transition {d = 1, z = 1), this scaling yields k ~ \fi — p, c \ 
which means that, although k ^ as /) -> l/2~,itis not 
a power law. 

In summary, we have presented QMC results showing 
that, contrary to previous results [llj where U > 2V, 
the extended one-dimensional soft core bosonic Hubbard 
model does have SS phases when doped above half filling 
when both U and V are large such that the system fa- 
vors multiple rather than near nighbor occupancy. From 
the scaling of f2(fc), p s and k, we calculated the dynamic 
critical exponent, z, for all the transitions. We showed 
different universality classes due to the different behavior 
of fl(k) depending on whether the CDW is approached 
from below or above half filling and we highlighted sim- 
ilarities and differences with the two dimensional case. 
Finally we note that lo! SS phases have been predicted 
for commensurate mixtures of two bosonic species, each 
near the hardcore limit, but where the interaction be- 
tween the species is not too large |18|. 
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